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ABSTRACT 


Velocity fluctuations over ocean waves and wave heights 
are examined using joint probability density function-condi- 
tional mean function (JPDF-CMF) analyses. Results are conm- 
pared with theoretical predictions from wind-wave coupling 
models which provide for interaction between the wave in- 
duced motion and turbulence in the airflow, and also with 
predictions from potential flow theory. 

Phase and amplitude results for wave-related fluctua- 
tions are examined in detail and consistent departures from 
potential flow are noted. These departures are phase differ- 
ences in between u and w which are not in quadrature, ampli- 
tudes of u and w which vary throughout erroneous’. and 
phase shifts of u and w which differ from those predicted 
by potential flow theory. 

Conclusions for the role of turbulence in the interac- 


tions are based on JPDF-CMF and phase amplitude analyses. 
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IL. INTRODUCTION 


More detailed studies of the near atmospheric boundary 
layer structure and energy exchanges are necessary if solu- 
tions to several applied geophysical problems are to be 
solved. Some of these applied geophysical problems include: 
the forecasting of surface wave conditions which require 
wind speeds at specified layers in the surface layer for 
semi-empirical formulae, and specifications of kinetic 
energy losses and thermal energy gains in the surface layer 
that are necessary for extended period weather predictions. 

Previous methods for obtaining boundary fluxes and deter- 
Mining the structure of the air layer near the ocean surface 
consisted primarily of applying empirical relationships, 
which have been found to be valid over land. These previous 
efforts generally made use of the assumption that interac- 
ELons associated with wind-wave coupling are negihigibhe. 

"An example of the need for a better understanding of the 
atmospheric boundary layer structure was the discrepancies 
in the interpretations of results presented at the Confer- 
ence on the Interaction of the Sea and Atmosphere at Fort 
Lauderdale, Florida in December 1971. Various results were 
marked by conflicting interpretations when examined with 
Strictly empirical approaches. Most evident in these were 
the various representations and values obtained for the drag 


coefficients as well as differences from intercomparisons 





between stress estimates obtained from different methods. 
The basic assumptions in obtaining these results were the 
form of the wind profile and neglect of significant depend- 
ence on the waves. General descriptions of the influence 
of wind-wave coupling within the framework of theoretical 
models combined with observational data should yield re- 
sults which delineate reasons for these discrepancies. 

It is appropriate to perform Vee oe analysis of velo- 
erty eon over natural waves at this time not only 
because of the need to consider the dynamics of wind-wave 
coupling to explain discrepancies among results but also 
because of the present stages in both theoretical and obser- 
waeromeal efforts im thismareay, Recent adwances inwetieoreti- 
cal and observational studies have brought both to the 
point where the nesults should and can be redated. 

Theoretical studies on wind-wave coupling have pro- 
gressed from the initial quasi-laminar descriptions (see 
Miles 1957) wherein the wave-induced motion interacts only 
with the mean wind to those wherein interaction occurs with 
the turbulence in the airflow (Yefimov 1970). Solutions 
from theoretical models which allow for turbulence in the 
airflow have revealed several features recently observed in 
results obtained over natural waves. These features in- 
clude amplitude and phase variation of the wave-related 
fluctuations with height and contributions by the wave- 
induced motion to the cospectra of u and w. Similarities 


in these results suggest that theoretical models could be 





refined and the observational results interpreted better if 
the two types of results were intercompared. 

In this study, velocity fluctuations over natural waves 
will be examined in order to show that the airflow does not 
necessarily follow the potential flow definition inherent in 
early wind-wave coupling models. For this, recent theories 
and discussions of wind-wave coupling will also be con- 
Sidered. These include those by Yefimov (1970), Manton 
(1971), Davis (1972), and Davidson (1970). 

The data examined represent simultaneous Senn mea- 
Smements at different lewels over oceam waves on the Atilan- 
tic Ocean. The measurements were made from the Scripps 
Institution of Oceanography Floating Instrument Platform 
(FLIP) which appears in Figure 1. The measurements were 
made during the Barbados Oceanographic and Meteorological 
Experiment (BOMEX). The BOMEX experiment was a joint inves- 
tigation of several agencies and universities for the pur- 
pose of establishing the energy budget of an atmospheric 
volume 500 km on a side and about 5 km high. The volume was 
located in the Trade Wind region of the Atlantic Ocean as 
shown in Figure 2, The data represent measurements near the 
center of the volume, in the vicinity of 16°N latitude and 
57° longitude, 

In this study, joint probability analysis on velocity 
and wave are the primary source for information on inter- 


actions that could have taken place due to wind-wave 
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Figure 1. Side and Top View of FLIP. 
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Figure 2. BOMEX and FLIP Location During 19-28 May 1969. 
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coupling. From the analyses, interpretations will be made 
on phase-amplitude relations, non-linear properties and 
Panerimerenron Ore tie ™type of tirbublencen™present. It iS siwg- 
gested that these statistical procedures will find wide- 
spread acceptance in the study of micro-scale meteorological 
PPeismetia. 

Data examined in this study were also analyzed by spec- 
tral analyses in a parallel investigation by Bingham (1972). 
Some of his results will be referenced during the discus- 
Sions. Spectral analysis, however, is limited when applied 
to studies involving phase-amplitude relationships between 


wind and wave. 
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Il. ASPECTS OF THEORETICAL AND OBSERVATIONAL RESULTS 


A primatyemotivatwen for the design ofethis investiga- 
tion was the close agreement between recent observational 
results and numerical solutions from recent wind-wave coupl- 
ing models by Yefimov (1970) and Reynolds (1968). In for- 
mulating their models, Yefimov and Reynolds allowed for 
interaction between the wave-induced motion and the turbu- 
lence in overlying shear flow. Although most theoretical 
formulations of the linear theory for wind-wave coupling 
have several common assumptions; such as, two dimensional 
Flow and a specified mean wind dependent only on height, 
and common simplifications; such as a single component wave 
faledkd and suitable linearized forms for the defining equa- 
tions, a distinguishing feature of recent formulations has 
been the role prescribed for turbulence in the airflow. In 
Conmitmast, Mishes .(1957), in Pe icikat Mie an early quasi- 
laminar model, neglected turbulence except for its role in 
determining an assumed logarithmic mean profile. 

Davis (1972) has shown that anistropic stresses, which 
amewmost. likelyeto bemsturbulent Reynolds stresses, are 
essential to the process of energy flow to the wave. These 
stresses, however, have been difficult to isolate in the 
various experimental tests of the models. Careful measure- 
ments and improved analyses of the fluctuating velocity 
components and the turbulent stresses are apparently what 


are required. 
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A. DISCUSSION OF WIND-WAVE COUPLING THEORY PREDICTIONS 

Miles' theory is a quasi-laminar wave generation mechan- 
ism in which the wave-related motion is assumed to follow 
the potential flow theory. The key to Miles' formulation 
is the existence of the critical level which corresponds to 
the level where the wind speed is equal to the wave speed. 
A schematic of Miles' assumed airflow is shown in Figure 3. 
The wave in this figure is stationary. Below the critical 
level Gigeceamiaimen the air appears to move to the left 
since it moves at a slower speed than the waves. Above the 
critical level, the wind moves to the right. Miles' theory 
is based on the premise that waves only modulate the mean 
mea Clow over them and that turbulence is™not™ influenced. 
Reynolds stresses are then dependent only on the height 
above the mean wave surface. Phase and amplitude features 
of the two components of velocity, and the stress uw, with 
respect to the wave, as predicted by potential flow theory 
are illustrated in Figure 4. 

Recently, Yefimov (1970) and others have allowed for 
turbulence in the airflow, which requires taking into 
account nonlinear terms in the equations of motion. A way 
to approach this has been to define the total vertical 


transfer of horizontal momentum as, (Yefimov 1970), 


du Ju 
u'w k : + ie as 


ow 


yz) 


where u' and w' represent total velocity fluctuations, u and 


w are Wave-related fluctuations, u is the mean wind, kA = 
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Figure 3. Streamlines of Sinusoidal Perturbations 
in a‘Shear Flow Near the Critical Level. 
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Figure 4. 


Phase-Amplitude Relations Predicted by Potential 
Flow Theory for the Two Components of Velocity 
(u and w) and the Stress uw, with Respect to the 
Wave. 
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mime snd 7 is anwexchange c@efiicient™melated to the wave- 
induced motion. 

Several features in Yefimov's solutions are due solely 
to the prescribed interaction between the wave-induced 
motion and turbulence. As such, they occur independent of 
the height or even existence of a critical level. In conm- 
parison, all nonpotential features of the airflow in initial 
quasi-laminar formulations were intimately connected to the 
dynamics at arte "eritical level'. These features were thie 
variations with height of the phase difference between the 
velocity components and the wave. Phillips (1966) and 
Kinsman (1965) describe the predicted phase shift for the 
Weweeical velocity component as being from over the downwind 
node toward the crest as the ‘critical level' is approached 
from the surface. Stewart (1969) and Reynolds (1968) were 
able to verify results from a wind-wave tunnel experiment 
With phase predictions from a nonlinear formulation for 
Wind-wave coupling and their work will be referenced for 
this aspect. 

Consideration of results with respect to theoretical 
approaches for the consequences of wind-wave coupling will 
be confined to prescribed or predicted characteristics of 
the velocity field in these models. These characteristics 
include: 

a. Variation of the amplitude of wave-related fluctua- 


tions with height; 
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pe. Contributions (negative or positive) to the 
momentum flux; and 

c. Variation of the phase difference between wave- 
related ex icjesiecaiacrien mao and the wave with 


height. 


B. STATISTICAL CONSIDERATIONS AND PROCEDURES 

Both joint probability density function-conditional mean 
function analyses (JPDF-CMF) and spectral analyses will be 
used to identify properties of fluctuation velocities asso- 
ciated with the waves' presence but only the JPDF-CMF analy- 
ses will be described in some detail. 

The following diagram outlines wwaseduaaa— used in the 
JPDF-CMF analyses. The procedures will be discussed in the 


ensuing paragraphs. 









Normalization and Computation 
Of Preliminary Stattstics of 
the Variable 





Y 


| Band Pass Filtering | 


Y 


| Computation of JPDF-CMF | 


+ 


Computation of Phase-Amplitude 
Information 


ilies PrRePiminary Statistics Computation and Normalization 


The procedures were applied to normalized values of 
the original data. Normalization was accomplished. by 


mii statistics were compiled and plotted at the W. R. 


Church Computer Facility of the Naval Postgraduate 
School utilizing the IBM computer and a CALCOMP plotting 
device. _— 
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dividing the deviation from the mean by the standard devia- 
tion; for example (x - x)/o. 

Z. Bande Pass Filtering 

A band pass filter was applied in order to focus 

the analyses on fluctuations near the frequency band corres- 
ponding to predominant surface waves. Specifically, 
application of the band pass filter allowed analyses of 
wave-induced fluctuations to be clear of any high or low 
frequency contamination. A numerical inverse transform 
filter was applied and was designed on the basis of pre- 
vious spectral results. The response curve for the selected 
band pass filter appears in Figure 5, where low and high 
frequency cut-offs are observed to be approximately .1 and 


.8 Hz respectively. 


1.0 


J nm (Hz) ! 10 
Figure 5. Numerical Inverse Transform Filter. 
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3. JPDF-CMF Computations 

JPDF-CMF procedures adopted for this study have been 
fully described by Holland (1972). In general, a JPDF re- 
presents the probability of a joint occurrence for a pair 
of variables and the CMF represents the average value of a 
third variable in a joint class interval for all observa- 
tion times when JPDF variables had values within the joint 
interval. Figure 6 shows a Tae Ae nel array represent- 
ing the trivariate statistical relationship described by 
the JPDF and CMF. From this figure, the probability of 
joint occurrences in two variables can be compared with a 
corresponding conditional mean of a third variable. 

4. Phase-Amplitude Theory and Computation 

Available theoretical (Yefimov 1970, Manton 1971) 
and observational (Davidson ' 1970, Holland 1972, Davidson 
and Frank 1973) results indicate that one of the important 
features of interest in this type of study is the statisti- 
cal dependence of fluctuations in the airflow on the phase 
of the wave. Such relations were determined by Davidson 
and Frank (1973) using procedures which are used in this 
study. It is in regard to phase-amplitude information that 
greater latitude for interpretations can be made in JPDF- 
CMF analyses than those obtained by spectral analyses. For 
example, JPDF-CMF results should reveal nonlinear features 
Since the information is not a result of linear expansions 
of the fluctuations. Furthermore, the approach can be 


applied to comparisons between wave-amplitude classes. In 
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this study the features are examined with respect to small, 
large, and combined or total waves. 

Phase information was determined using a method des- 
cribed by Holland (1972) which is based on JPDF-CMF analyses 
procedures. Phase-amplitude information is possible by 
such procedures because polar coordinates which represent 
the phase angle and amplitude of a variable such as the 
wave height (n) can be determined from the JPDF array for 
the variable, n, and its time derivative, fh. As an example 
of this, u(n, n) would show the dependence of the horizon- 
tal component of wind velocity on the phase and amplitude 
of the waves. As will be seen in later discussions, the 
phase is defined in terms of octants measured counterclock- 
wise in the nh, n coordinate system (JPDF array) from the 
positive axis. JPDF-CMF arrays which contain phase- 
amplitude information for u,v,w,uv,uw, and vw with respect 
to the waves will appear in the discussions of a following 
section. The phase-amplitude results are easily summarized 
by octants for each period analyzed. 

In determining phase-amplitude fave onto a para- 
meter, TV; called the ‘statistical eddy period' is also of 
interest. To is defined by the ratio 


Oo 
tT = 25 — 
oO 


Os 


where o and 6 are standard deviations of the variable and 
its derivative. If the statistical eddy periods compared 


differ greatly between the independent and reference 


23 





variable, phase-amplitude information for the reference 
variable has a greater possibility of only representing the 
mean of numerous random fluctuations. In general, the 
statistical eddy periods agree as can be seen in Table I. 
With respect to possible influence by FLIP's motion 
on the results, analyses by Rudnick (1967) have indicated 
that observed fluctuations in u and w are 107 times larger 
than could be induced by FLIP'’s Bee sige. Therefore, FLIP’s 


motion will not enter into the discussion of the results. 
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Tit. DESCRIP TLON OF MEASUREMENTS 


The data available for this study consisted of simul- 
taneous recordings of wind velocity components (with hot- 
film anemometers), air temperature fluctuations and wave 
heights. The levels of measurement ranged from 2 to 8 
meters above mean water level and the sensors and recording 
Systems were capable of measuring Fe lees and temperature 
fluctuations of frequencies up to 50 Hz. There were about 
54 hours of recordings in 40 separate observation periods 
to select from. On the basis of meteorological and wave 
conditions and on previous spectral analyses, four periods 
were chosen for the analyses. 

In the following paragraphs specific eee on the 
measurements and general conditions for every period will 


be described. 


A. FLOATING INSTRUMENT PLATFORM (FLIP) 

A discussion of FLIP has been provided by Bronson and 
Glosten (1968). The Floating Instrument Platform is a 
Super-stable, open-sea, free-floating platform, which is 
towed to a site while in a horizontal position and then 
partly flooded to cause the structure to flip to a vertical 
position for use as a research platform. Figure 1 included 
a side and top view of FLIP, while Figure 7 illustrates the 
working section of FLIP which remains above the water line 


when FLIP is in the vertical position. 
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Sensor Locations on FLIP. 
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In general, FLIP performed in a satisfactory manner 
throughout its use in BOMEX. Wind tunnel tests revealed 
that a minimum distortion to the natural flow of air past 
the intended location of an stiroara mounted vertical sen- 
sor supporting mast occurred when the deck of FLIP was 


positioned facing into the wind. 


B. EQUIPMENT AND PROCEDURES 
fr wensor Mounting Amrangement 

Figure 7 shows the locations of sensors for turbu- 
lence, profile and wave measurements on the vertical mast 
attached to FLIP. Sensors were mounted at the 2, 3, 4, 6, 
and 8 meter levels and all sensors were more than 1 meter 
windward from the vertical mast. As can be seen in Figure 
7, the vertical mast, supported by a horizontal boom, is 
about 15 meters crosswind from FLIP's hull. The mast. was 
attached to a carriage that could move horizontally on the 
boom. In addition, the mast could be turned on ‘ horpzon— 
tal axis so that it could be brought alongside a working 
platform to replace and adjust sensors during the experi- 
ment. 

No provision for adjusting the vertical position 
of the mast, independent of FLIP, were available. This 
made it almost impossible to maintain sensors in a mean 
level position. Provisions were made, however, to correct 
for this discrepancy. In general, it is felt that instru- 


ment mounting arrangements have not significantly affected 
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the accuracy of the data values with regard to the interpre- 
tations made. 
2. Wind Component Measurement 

Wind components were meaSured with hot-film, con- 
stant temperature anemometer oe t Hee with linearizers. 
Each Probes had three quartz-coated mutually perpendicular 
sensors 0.15 millimeter in diameter and 2.0 millimeters 
long. Each sensor is a glass rod covered with a platinum 
film onto which a quartz coating has been sputtered (see 
Figure 8). The sensors were electrically isolated so that 
simultaneous measurements could be made. The orthogonal 
array made it possible, with trigonometric relationships, 
to determine the instantaneous vector components of the 
wind from simultaneous measurements. The frequency response 
of these sensors, appeared to be well over 50 Hertz in the 
conditions under which they were used. 

The constant temperature anemometer system measures 
velocity by sensing the heat loss from the heated sensor 
to the moving air. If the air density is nearly constant 
and if the sensor is maintained at a constant and relative- 
ly high temperature, the heat loss is dependent almost 
entirely on air velocity. The sensor is part of a bridge 
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A Millimeter Scale 


Three-Sensor Hot Film Probe. 


Figure 8. 


Shown on the Left. 


is 


3:0 





3. Wave Measurement 
Wave data were obtained from a resistance gauge 
Made available by Dr. R. E. Davis, Scripps Institution of 
Oceanography. The probe consisted of a non-conducting 
tube, approximately 2.5 centimeters in diameter, with a 
conducting wire wrapped spirally around it. It was posi- 
tioned about 5 meters inboard from the vertical mast that 


held the velocity and temperature sensors. 
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IV. PRESENTATION OF RESULTS 


ite Tresurts are irom four different periods. The per= 
iods were in three separate days in May 1969: 19 May (1613 
-1735 GMT), 27 May (0311-0435 GMT), and 28 May (1756-1916 
and 2028-2150 GMT). Each period was approximately 80 min- 
utes in duration. Hereafter, these periods will be referred 
to as periods 1, 2, 3, and 4. Each period consisted of 
several subsets with regard to JPDF-CMF and phase-amplitude 
analyses; periods 1 and 2 consisted of four subsets of 20 
minutes each, and periods 3 and 4 consisted of three sub- 
sets of 20 minutes each. Separate reference will be made 
to individual subsets within these periods when required. 

Airflow data for periods 1, 2, and 4 were from the 8 
meter level and those for period 3 were from the 3 meter 
level. Data considered were the u, v, and w velocity con- 
ponents and wave heights. 

General conditions during the four periods are sum- 
marized in Table Il and Figure 9. General conditions 
during all periods could be described as being steady with 
respect to both waves and wind. Certainly the periods did 
not have active wave growth or rapid decay. These condi- 
tions are representative of conditions over most of the 
oceans when swell propagates into regions of relatively 
light winds. Therefore, these results are relevant to des- 


cribing the surface layer over the ocean. 
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Spectral results from these periods have been described 
by Bingham (1972). They show corresponding extrema in the 
velocity variance, covariance spectra and wave spectra and 
reflect a few characteristics of the predicted solutions, 
as will be noted in the results. UW cospectra had extrema 
which, in fact, represented either enhanced or decreased 
momentum transport when compared to that at frequencies not 
associated with wave spectra extrema. Potential flow does 
not predict such extrema in the uw cospectra because in 
theory, u and w are in quadrature. Such contributions, 
both negative and positive, are interpreted as being those 
predicted by the —r theoretical solutions. 

Physical processes can be realized from the arrays re- 
presenting probabilities of joint occurrences of fluctua- 
tions of various magnitudes in (u, w, uw, etc.) in 
conjunction with an average surface elevation during the 
joint occurrences. 

JPDF-CMF patterns for potential flow are presented in 
Figure 10. As indicated in Figure 4, potential flow 
theory predicts a maximum positive value of uw occurring 
in quadrant II and IV. The JPDF-CMF pattern in Figure 10, 
(Panel A), represents such a relationship. Similarly, u 
and w relationships which appear in Figure 4 are repre- 
Semted in Panels B and F. Pamel C and E will be discussed 
in the following paragraph. Panel D reflects the absence 


of any correlation between the CMF, v, and the JPDF (yj, n). 
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Panel C represents the probability of joint occurrences 
of fluctuations of various magnitudes of u and w with res- 
pect to the wave height. The following diagram interprets 
the four quadrants in Panel C with respect to momentum 
transfer. Panel E has, of course, the CMF, wave derivative, 


rotated 90 degrees clockwise from its orientation in Panel 





Cc. 

+ W 
Upward transfer of Upward transfer of 
deficit momentum excess momentum 
ij ar U 
Downward transfer of Downward transfer of 
deficit momentum excess momentum 

- W 


me fol herime seetions imetude detailed discussions of 
features in JPDF-CMF patterns for each of the four periods 
commencing with period 1 and continuing through period 4. 


A summary will follow the presentation of results. 


i ~fHERIOD 1 

Significant features identifiable within JPDF-CMF re- 
sults from this period are evidence of mechanical turbulence 
associated with hydrostatic stability, decreased momentum 
transfer due to the waves' presence, large deviations from 
potential flow and apparent changes in u and w phase rela- 
tions with respect to the wave due to height variations of 


the critical level. 
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With regard to the evidence of mechanical turbulence, 
the JPDF (u,w) patterns can be interpreted for these fea- 
tures in the fluctuations. Holland (1968) proposed that if 
turbulence is primarily mechanical and is occurring under 
stable conditions then a wider distribution of both u and 
w should occur at positive u than at negative u and the 
main mode, at high probability should be shifted towards 
negative u. Distributions which fit such a description 
appear in the JPDF pattern in Figure 1l, Panel C. If the 
hydrostatic conditions had been unstable, the opposite 
effect, a narrow distribution of u and w at positive u and 
the main mode shifted towards positive u would have been 
seen. 

JPDF patterns below the 5 percentage level in Figure ll, 
(Panel C) are slightly elliptic, asymmetrical with a major 
axis 135 degrees counterclockwise from the positive u-~axis. 
Statistically, a small probability of joint occurrences of 
fluctuations of various magnitude of u and w is observed 
in quadrants II and IV, indicating that, due to the waves' 
presence, upward transfer of deficit momentum and downward 
transfer of excess momentum occurred. 

From cospectra in Figure 12, it is observed that upward 
transfer of momentum was evident near the wave spectra 
peak. Similarly, the center patterns in Figure 1l (Panel 
C) are nearly circular and correspond to decreased downward 
transfer as shown in the cospectra, Figure 12. Recall that 


in these interpretations the band-pass filter (see Figure 5) 
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included frequencies of a greater range than the wave 
spectra peak band. 

With respect to phase relations, large deviations from 
Botential flow can be noted by comparing u(n, n) and w(n, 
N) patterns of Figure 13, (Panels B and F) with comparable 
panels in Figure 10. In particular, a rotation of about 
45 degrees is seen in the CMF of u(y, nN). This is an 
obvious shift from that predicted hgeeporen erie flow (Figure 
10, Panel B). CMF patterns in w(, n) are not as dominant 
as the patterns to be considered later for period 3, yet a 
rotation of about 45 degrees can be discerned. The phase- 
amplitude relationship in Figure 14 further reflects the 
above deviation from potential flow. The phase-amplitude 
representation shows w to be nearly in-phase with the u 
component. In contrast, potential flow theory predicts 
that u and w will be in quadrature as was shown in Figure 4. 

Karaki and Hsu (1968) observed phase results similar to 
those just described by comparing u and w perturbation over 
water waves in a laboratory experiment. These were ob- 
served, from their analysis, near the region of the 
critical level and were due to a shift in the phase of w 
relative to both the wave phase and u~phase. The relative 
u-w phase was observed to be Near zero an eee critical 
level. Stewart (1970) also described results where the 
phase shift of both u and w with height resulted in u and w 


being out of quadrature and nearly in-phase. 
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I-DE=aCite Pathernms for Periodel, Subset 3. 
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Figure 14. Phase-Amplitude Results (Components) for 
Period 1, Subset 2. 
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The relations just described can be compared to those 
Securring later in the period by comparing the results in 
Figure 14 to those in Figure 15. Results in Figure 15 cor- 
respond to a time interval which was 42 minutes after the 
interval corresponding to the results in Figure 14. 

Considering Figure 15, the phase difference between u 
and w gradually increased as the period progressed until a 
shift of nearly 45 degrees was reached during the last sub- 


set of the period, Figure 15. 


Be PERIOD 2 

Significant features observed in this period were: 
evidence of mechanical turbulence associated with hydro- 
static stability, apparent potential flow conditions and 
variations in the magnitude and direction of the wave-in- 
duced momentum transfer. Also changes appear to have 
occurred during the period in the physical make-up of the 
wave. This period, unlike any other period, had near 
mecential flow conditions in both u and w. 

Pithmrespec: Co the rohegof mechanical turbulence, im 
Figure 16 (Panel C) n(u,w) JPDF patterns are highly ellip- 
tic with asymmetrical rings orientated so that the major 
axis is 135 degrees from the positive u-axis. The highly 
elliptic pattern will not be observed in any of the other 
periods. Again the turbulence is interpreted to be 
mechanical as in period 1 on the basis of the shift of the 
JPDF pattern. The corresponding cospectra in Figure l1/ 


show only slight enhanced momentum transfer. 
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Phase-Amplitude Results (Components) for 
Period 1, SubSet 4. 
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Figure 16. JPDF-CMF Patterns for Period 2, Subset 4. 
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Mmemre 1/7. opectral Results (Period 2) Including Co- 
Vanlance Specira and’ Wave opece ia. 
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With respect to similarities with potential flow pre- 
dictions, u(n,n) and w(n,n) results in Figure 16 (Panels 
B and F) and throughout the period reveal no significant 
shift of the CMF patterns from those predicted for poten- 
tial flow (Figure 10). Also, w(n,n) CMF patterns are weak- 
ly defined when compared to other periods. This indicated 
a smaller probability of joint occurrences in average 
fluctuations of :. With Smead 1. 

With respect to the suggested possibility of a differ- 
ent wave shape, JPDF patterns in Figure 16 (Panels B, D, 
and F) have a marked difference from those in any other 
period, such as those in Figure 11. Where previous and 
later periods have nearly concentric patterns showing only 
slight deviations from those of a sinusoidal wave, this 
period has more-elliptical patterns with a slightly in- 
creased gradient over the wave crests. The corresponding 
Metatistical eddy period’ is also much lower during this 
period than the other periods (Table I). This could be 
indicative of a marked change in the physical make-up of 
the wave, waves having sharper crests and broader troughs. 

With respect to momentum transfer, uw, in relation to 
amplitude classes, the large wave results in Figures 18, 
H9, 20, and 21 appear to shift from an upward transfer of 
deficit momentum to a downward transfer of deficit momentum 
as the period progresses. To a ee degree, the shift 
is also present in the total wave. This is interesting in 


view of the results in Figure 4 which show the general wind 
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Figure 18. Phase-Amplitude Results (Stresses) for 
BeLTlod 2. oupaet 12 
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Figure 19. Phase-Amplitude Results (Stresses) for 
Period 2, Subset 2. 
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Figure 20. Phase-Amplitude Results (Stresses) for 
Peamlode 2, Swbiget: 3m 
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Figure 21. Phase-Amplitude Results (Stresses) for 
Period 2, Subset #4. 
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conditions for all periods. The wind has been increasing 
with a large increase during the last part of this period. 
Melis can possibly be attributed to what Yefimov (1970) 
called the complex phase as the wind speed approaches the 
phase speed of the predominant surface wave. 

The last feature noted for this period is that the wave- 
form of the wave related stress (uw), Figures 18, 19, 20 
and 21, also varies considerably over the period. Consider- 
ing the large wave amplitude class, early in the period the 
largest amplitude exists in quadrant II but as the period 
progresses, the dominent quadrant shifts to quadrant I. 
With respect to the total waves, however, initially the 
wave-form is similar to the potential flow prediction but 
quadrant II becomesdominant later in the period. With res- 
pect to the third subset, the large wave amplitude class 
has the maximum value in quadrant III as does the total 
wave amplitude class. The final minutes of this ‘period 
continue with the same relation found in subset three, that 


of maximum amplitudes in quadrant III. 


Ge PERIOD 3 

This period differed from the other periods in that 
measurements were made at 3 meters. Significant features 
in period 3 appear to be reduction in excess momentum 
transfer at the frequency of the wave spectrum peak, dis- 
tinct deviations from potential flow in the phase of w, and 
uw phase-amplitudes which differ significantly from potential 


flow theory. 
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With respect to momentum transfer, JPDF patterns in 
Figure 22, (Panel C) have the appearance of concentric cir- 
cles. They do have a slight asymmetry in quadrant IV, re- 
presenting some downward transfer of excess momentun. 
Spectral results for this period appear in Figure 23. The 
uw cospectral results reflect excess momentum transfer at 
the frequency of the wave spectrum peak. Such a spectral 
result does not agree with the JPDF axis rotation in Panel 
C, which is approximately 315 degrees counterclockwise from 
the positive u-axis. 

The CMF patterns are now examined in order to determine 
the role of the waves in producing the (u,w) JPDF patterns 
in Figures 22 and 24 (Panel C). For potential flow the 
CMF patterns for the waves would have maximum positive 
values associated with negative u and zero w (Figure 10, 
Panel C). Maximum negative CMF values would be associated 
with positive u and zero w. A deus nee deviation from the 
pattern is the clockwise rotation of both maxima and minima 
in the CMF patterns of approximately 35 degrees in Figures 
22 and 24 (Pemel C). | 

With respect to phase relationships, the significance 
of the phase shift in both extrema can be evaluated further 
by examining the JPDF-CMF patterns for u(,n) and w(n,n) in 
Panels B and F. The CMF pattern orientation predicted by 
potential flow describes that for u fairly close with only 
a slight rotation as shown in Figure 22, Panel B. [In the 


case of w, however, where maxima/minima extremes are 
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Pigure 24. 


Ir -CMreratterns for Period 3, Subset 3. 
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supposed to be concentric around zero values of n.- There is 
a distinct rotation backward toward positive and negative 
values of 7. This also can be seen from the summarized 
phase-amplitude for this period, Figure 25. For all ampli- 
tude classes, u is also comparable to potential flow, but 
wis not. The latter is made evident by the minimum value 
of w occurring before the inflectien point on the wave. 

This represents, perhaps, atadence of downuwardmrnanster of 
excess momentum over the phase of the wave due to 
mechanical turbulence. 

To further emphasize the distortions that appear in 
Figure 25 for u and w, the uw phase-amplitude results and 
uw(, n) JPDF-CMF results will be considered. In potential 
flow, Figure 4, maximum values of uw occur behind the crest 
and following the trough, whereas minimum values occur 
immediately preceding the crest and the trough. The phase 
relations appearing in Figure 26 are comparable to those 
predicted; however, the amplitudes vary over the different 
phases of the wave. For example, quadrant IV amp Pi eude 
values are greater than the others. This is also apparent 
in the JPDF-CMF results. The latter result in some down- 
ward transfer of excess momentum. Quadrant Iii ampli- 
tudes are also greater than those in quadrant I which 
suggests that the occurrence of downward transfer, of 
deficit momentum predominates as associated with hydro- 
etatic stability, and, therefore, consistent resuits are 


obtained in this period. 
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Figure 26. Phase-Amplitude Results (Stresses)for 
Period ’3.- Subset. 1. 


60 





Considering the same features as above but for the JPDF- 
CMF uw(n, nn), this period has maximum values of uw, Figure 
22, (Panel A), occurring in quadrant IV but rotated slightly 
bammaras the trough position. This rotation, although not 
as apparent in the phase-amplitude results in Figure 26, is 
then consistent with the previously noted departures from 


potential flow theory. 


D. PERIOD 4 

The significant features noted in this last period are 
a small negative u-w correlation, a shift away from near 
potential flow for u and w, and definite stress deviations 
which suggsst wind-wave interactions. 

With respect to momentum transfer processes, n(u,w) 
Petterus im Migures 2/ and 28, (Panel C), appear as minor 
elliptic, asymmetric circles. This appearance signifies 
a smaller probability of joint occurrences of negative 
correlation of u and w. However, the inner rings in 
Figure 27 (Panel C), represent a slight elliptic pattern 
in the opposite direction and reveal a definite positive 
correlation. Similarly, the phase relations of u and w in 
Figure 29 show u and w to be nearly in-phase. In subset 
3 of this period, Figure 28, (Panel C), the positive corre- 
lation is not as evident suggesting a shift between u and 
w phases which can be verified by the phase-amplitude re- 


sults in Figure 30. 
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Figure 29. Phase-Amplitude Results (Components) for 
Period 4, Subset l. 
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Miemorlheieaeton Of CMF patterns in Figures 27 and 28, 
(Panel B), corresponding to u(y, n) are similar to those 
for near potential flow during the beginning of the period, 
then gradually rotate as the PNW roc 1 O86 s Eo a final 
rotation of about 30 degrees behind the crest. Similar 
rotations are not evident in the whem) “Gli patterns, 
(Panel F). Phase relationships in Figure 29, however, show 
maximum values of w occurring forward of the front node, 

a feature noted by Karaki (1968) and Stewart (1970). A 
possible reason for this rotation of the CMF pattern may 
be the turbulent interaction between the near surface wind 
and the wave. 

With respect to phase variations of the stress uw, 
initially the maximum positive values of uw, Figure 31 are 
forward of quadrants II and IV and are not equal in magni- 
tude. Later in Figure 33, the magnitude of uw is nearly 
equal in both quadrants II and IV. In both cases, however, 
a sufficient shift is noted from potential flow. This was 
also observed by Kendall (1970). The uw(n, n) patterns, 
Figures 27 and 28, (Panel A), also reflects these deviations 
from potential flow since higher CMF values appear in 
quadrants II and IV. However, looking at the phase rela- 
tionships in Figure 31, quadrant IV, downward transfer of 
excess momentum is the largest for both large and total 
waves. During the middle of the period, Figure 32, quad- 
rant II is predominant over large waves, while quadrant IV 


is the dominant over total waves. During the last part of 


65 









CM/SEC 





YCM/SEC 


Figure 30. Phase-Amplitude Results (Components) for 
Perieds4, Subset 3. 


66 








Ne] == 
va 
lax iN 
4 \\ \ 
UW i— : i | | = S, 
LF NZ 
M7 Nf 
f°’ J 
UV / \ / \ | as 
— _- a 4 a ] ra \ a 
—_—- 
 -~ . an 
n f— =, ———— 
ww f) \ ; 
ae / . 
\ / S 
eZ. 
A B C 
SMALL LARGE TOTAL 


Bieunge 31. 


Phase-Amplitude Results (Stresses) for 
Period 4, Subset 1, 


od 











lf oN 
7 Ne 7 
N 
F \ ‘ 
PAs 
/ \ 
Sf 7 ah Ss 
| \ 
y, 
“| 
B c 
SMALL LARGE TCTAL 
Figure 32. Phase-Amplitude Results (Stresses) for 


Period 4, Subset 2. 
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Figure 33. Phase-Amplitude Results (Stresses) for 
Pernod 4; SubDSet 3. 
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the period, Figure 33, large uw stress amplitudes appear in 
quadrant II for both large and total waves. In general, 
there is a shift from an initial downward transfer of 
excess momentum to an upward transfer of deficit momen- 
tum. This suggests that, due to the waves' influence, 
momentum is being transfered towards the surface during the 
first half of the period and transfered away from the sur- 
face during the last half of the period. During the 
reversal of Heenan transfer between the first part of 

the period to the last part of the period, mechanical tur- 
bulence may well be the important factor in the momentum 


transfer. 





V. SUMMARY AND CONCLUSIONS 


Analyses were performed on observational data obtained 
over ocean waves. These analyses provided results on the 
phase and amplitude relationship of the wave-induced fluc- 
tuations. These results were interpreted with respect to 
potential flow predictions for a wave propagating beneath a 
shear flow and the role of turbulence in wind-wave coupling. 

A significant feature was the consistent rotation of 
Eve Clee pattern in the u(n, n) and wn, n) results from the 
orientation which would correspond to potential flow. This 
rotation was noticeable in all periods analyzed to some de- 
gree. Also observed was the absence of the quadrature rela- 
pLonscnip or u and w Ae predicted by potential flow theory. 

In comparison to potential flow predictions which have 
equal magnitudes of uw phase-amplitude in all quadrants, 
the uw phase-amplitude results revealed variations of magni- 
tude over different phases of the wave. This suggests that 
nonlinear processes are responsible for momentum transfer 
in wind-wave coupling. 

Therefore, wave~caused motion of the air in the near- 
water layer was observed to be essentially nonpotential. 
This was most evident in periods 1, 3, and 4 with a lesser 
deviation from potential aes during period 2. However, 
during period 2, a marked change was taking place in the 


physical make-up of the wave. 
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